Phase Behavior of Type-II Superconductors with Quenched Point 
Pinning Disorder: A Phenomenological Proposal 



Gautam I. Menon* 
The Institute of Mathematical Sciences, C.I.T. Campus, 
Taramani, Chennai 600 113, India 
(February 1, 2008) 

Abstract 

A general phenomenology for phase behaviour in the mixed phase of type- 
II superconductors with weak point pinning disorder is outlined. We propose 
that the "Bragg glass" phase generically transforms via two separate ther- 
modynamic phase transitions into a disordered liquid on increasing the tem- 
perature. The first transition is into a glassy phase, topologically disordered 
at the largest length scales; current evidence suggests that it lacks the long- 
ranged phase correlations expected of a "vortex glass". This phase has a 
significant degree of short-ranged translational order, unlike the disordered 
liquid, but no quasi-long range order, in contrast to the Bragg glass. This 
glassy phase, which we call a "multi-domain glass", is confined to a narrow 
sliver at intermediate fields, but broadens out both for much larger and much 
smaller field values. The multi-domain glass may be a "hexatic glass"; alter- 
natively, its glassy properties may originate in the replica symmetry break- 
ing envisaged in recent theories of the structural glass transition. Estimates 
for translational correlation lengths in the multi-domain glass indicate that 
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they can be far larger than the interline spacing for weak disorder, suggest- 
ing a plausible mechanism by which signals of a two-step transition can be 
obscured. Calculations of the Bragg glass-multi-domain glass and the multi- 
domain glass-disordered liquid phase boundaries are presented and compared 
to experimental data. We argue that these proposals provide a unified pic- 
ture of the available experimental data on both high-Tc and low-Tc materials, 
simulations and current theoretical understanding. 
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I. INTRODUCTION 



In the mean-field phase diagram of a pure type-II superconductor, Meissner, mixed and 
normal phases are separated by continuous phase transitions associated with lower (Hci(T)) 
and upper (Hc2(T)) critical fields |l|]. In the mixed phase, an applied magnetic field H 
enters the sample in the form of singly quantized lines of magnetic flux. These lines repel 
each other, stabilizing a triangular line lattice. Thermal fluctuations melt this lattice via 
a first-order melting transition, subdividing the domain of the mixed phase into solid and 
fluid phases of flux lines 0-0. 

With H = Hz, flux-lines in the lattice phase are parametrized by two-dimensional coor- 
dinates Ri, where Rj describes the position of the i*'^ line. Deviations from this arrangement 
define a two-component displacement field u(Rj, z), with an associated coarse-grained elastic 
free energy cost 

Fei = J dv^dz [^(Vx ■ u)2 + ^(Vx X u)2 + ^(9,u)^]. (1.1) 

Here u{v±,z) is the displacement field at location (r^,z). This term represents the elastic 
cost of distortions from the crystalline state, governed by the values of the elastic constants 
for shear (cee), bulk (cn) and tilt (C44). 

Quenched random pinning destabilizes the long-ranged translational order of the flux- 
lattice phase Such disorder adds 

Fd = J dr±dz Vrf(r_L, z)p{r±, z), (1.2) 

to Fei above, where represents the interaction of vortex lines with a quenched disorder 
potential Vd and the density field p(r_L, z) is defined by 

p{r^,z) = Y,5'-^\r^ - R, - u(R„^)). (1.3) 

i 

If the disorder derives from a high density of weak, random point pinning sites, it is 
convenient to assume that Vd is drawn from a Gaussian distribution and is correlated over 
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a length scale ^, the coherence length i.e. [Vd{x)Vd{x')] = K{x — x'), with K{x — x') a 
short-ranged function of range ^ 0. The notation x denotes {r±,z). 

Much theoretical attention has been devoted recently to the understanding of the sta- 
tistical mechanics defined by Eqns. and ( |1.2| ). The problem is the computation of 
correlation functions such as 

i?(rx,^) = [(u(rx,z)-u(0,0))2], (1.4) 

where the brackets (■) and [■] denote thermal and disorder averages respectively. At least 
two novel glassy phases arise as a consequence of such disorder In the relatively more 

ordered of these phases, the Bragg glass (BrG) phase, correlations of the order parameter 
for translational correlations decay as power laws |f^-[T2|: 



Cg(xx) = (exp^[G ■ (u(x^,z) - u(0,^))]) ~ l/lx^^^, (1.5) 



where G is a reciprocal lattice vector and f/c a non-universal exponent 1TT||T31 • "^^^ anoma- 
lously large translational correlation lengths inferred from magnetic decoration experiments 
on BSCCO [O and the resolution limited Bragg peaks obtained in neutron scattering from 



the mixed phase in this material at small H, support the existence of a Bragg glass 
phase in type-II superconductors. Hall probe measurements [l^ see a single first-order. 



temperature-driven melting transition out of the Bragg glass at such fields [p!^-[T7|l. 

At increased levels of pinning, the Bragg glass is believed to be unstable to a disordered 
phase in which translational correlations decay exponentially at sufficiently large distances 
PJTOd I.e. lima;^oo Cg(x) ~ exp(— /i|x^|). Such a phase must be separated by a phase 
transition from the Brass glass phase. At large H, experiments |jl8[ indicate a continuous 



transition from an equilibrium fluid phase into a highly disordered glassy state when the 
temperature T is reduced from the mean-field T(.{H) p^ , |20| . This disordered low-T, large- 
H state may be a new thermodynamic phase, the "vortex glass" (VG) phase, separated 
from the equilibrium disordered liquid (DL) by a line of continuous phase transitions [Q. 
However, the experiments may also be indicating a non-equilibrium regime A "frozen 



liquid" similar to a structural glass and an "entangled liquid" analogous to a polymeric glass 
have been suggested as alternatives to the equilibrium vortex glass pT| , p2| . 



The phase diagram of Fig. 1 summarizes the currently popular view of the phase be- 
haviour in the mixed phase of weakly disordered single crystals of anisotropic high-Tc systems 
such as BSCCO P,p!^, P^ , pil| ,p3|-p5[ . (This disorder is assumed to arise from a large density 
of weak point pinning sites in this paper |2^.) Fig. 1 implies that the BrG phase is unstable 
to the VG phase on increasing H. This feature follows from the belief that increasing H at 
a fixed level of microscopic disorder is equivalent to increasing the effective level of disorder 



27| . Fig. 1 also shows the Bragg glass melting directly into the DL phase in the intermedi- 



ate field regime (solid line). The VG to DL transition begins where the first-order BrG-DL 



transition line terminates; this transition has been suggested to be continuous [^. Recent 



m- 



Josephson plasma resonance experiments and magneto-optic studies in BSGCO p8|j29 
dicate that the underlying field-driven BrG-VG transition may actually be a discontinuous 
one . 

The importance of thermal fluctuations in the cuprates ensures that disordered phases 
such as the DL phase occupy much of the phase diagram in high-Tc materials In 
contrast, for most low-T^ systems, such disordered phases occupy a relatively small regime 
just below Tc(-ff). The phase behaviour in the mixed phase of low-Tc superconductors has 



been the focus of recent attention [pl| -|35| . These studies suggest a universal phase diagram 



for weakly pinned low-T^ type-II materials [^. Interestingly, this phase diagram appears 
to differ significantly from Fig. 1. 

The low-Tc systems studied have relatively large values of the Ginzburg number, 

= {kBTjHleef/'^. (1-6) 

which measures the relative importance of thermal fluctuations Here ^ is the coherence 
length, e the mass anisotropy, Tc the superconducting transition temperature and the 
thermodynamic critical field. In most low-Tc materials Gi ~ 10"^, whereas Gi ~ 10^^ for 
high-Tc superconductors such as YBa2Cu307_5 (YBCO) and Bi2Sr2CaCu208+<5 (BSCCO) 



For the layered dichalcogenide 2H-NbSe2 (T^ — 7.2K), anisotropy, weak layering, small 
coherence lengths and large penetration depths enhance the effects of thermal fluctuations, 
yielding Gj ~ 10^"^ . For the C15 Laves-phase superconductor CeRu2 (Tc — 6. IK), 



Gj ~ 10 ^ ||5S|, while for the ternary rare-earth stannides Ca3Rh4Sni3, (T^ — 8K) and 
Yb3Rh4Sni3, (T^ ~ 7.6K) , G^ ~ 10"^ For the quarternary borocarbide YNisBaC 



(Tc ~ 15. 5K), Gi ~ 10" |^2[- These systems can be made relatively pure: the ratio 
of the transport critical current jc to the depairing current jdp {jdp = cHc/SVQttX, with 
He = ^0/2^/2,11 \^) in these low-Tc materials can be as small as 10~^ or 10~^. For the 
high-Tc case, in contrast, jc/jdp ~ 10"^. These values of Gi (1-4 orders of magnitude larger 
than for conventional low-T^ materials) and of jc/jdp imply that the phase behaviour of 
weakly pinned, low-T^, type-II superconductors can be studied using such compounds, in 
regimes where the phases and phase transitions discussed in the context of the cuprates are 
experimentally accessible. 

Fig. 2 presents the proposal of this paper for a schematic phase diagram for type-II 
superconductors with point pinning disorder. This phase diagram incorporates and extends 
recent suggestions for the phase behaviour of low-T^ systems, proposed by this author and 
collaborators in Ref. It is argued that Fig. 2 should be a universal phase diagram for 
all type-II superconductors with quenched point pinning disorder, both low and high-T^, in 
a sense made precise in this paper. 

In Fig. 2, the term "multi-domain glass" (MG) describes the sliver of intermediate glassy 
phase, which we propose should generically intervene between BrG and DL phases. This 
phase is argued to be an equilibrium phase. At present, neither experiments nor simulations 
provide a clear indication of the nature of this phase nor constrain the various proposals 
made in this regard. One possibility is that the symmetry breaking which occurs across 
the MG-DL transition line is related to the replica symmetry breaking proposed recently in 
the context of the structural glass transition [^. A hexatic glass is an alternative which is 
attractive in many respects p4 |. 



Modern replica-based approaches |^ to the glass transition in the context of fluids with- 
out quenched disorder, which argue for a one-step replica symmetry breaking transition out 
of the disordered liquid phase, have recently been generalized to study the phase diagram 
of hard spheres in the presence of a quenched random potential |^5[ . Such a phase diagram 
bears a fairly close resemblance to phenomenological phase diagrams of disordered vortex 
systems, including the existence of both flrst-order and continuous transitions out of the 
disordered liquid |^5|. The nature of symmetry breaking here is subtle and it is unclear 



whether experiments would be able to distinguish between a genuine equilibrium glass tran- 
sition of this type and a non-equilibrium "freezing-out" of structural degrees of freedom. 
Nevertheless, in principle, signals of such an equilibrium transition should be visible in ther- 
modynamic measurements. Should these ideas be applicable to the disordered mixed phase, 
a discontinuous freezing of the disordered density conflguration characterizing the glassy 
phase would be expected, as against relatively smooth behaviour of the entropy and internal 
energy across this transition. 

It is believed that the lower critical dimension for phase coherence of the type assumed 
in models for vortex glasses exceeds three EH]; this would rule out the phase- coherent vor- 



tex glass 10 as an alternative. However, the absence of such long-ranged superconducting 
phase coherence does not rule out the possibility of various levels of order and associated 
symmetry breaking transitions within the vortex line assembly. It is to be emphasized that 
the precise pattern of this symmetry breaking is irrelevant to the phenomenology embodied 
in Fig. 2, in particular to the proposal of a generic two-step transition out of the BrG 
phase which occurs via an intermediate equilibrium phase with glassy attributes, in which 
length scales for translational correlations in the flux line array can far exceed those in the 
equilibrium fluid. Even if the true ground state were an equilibrium vortex glass in the 
sense of the original proposal, the topology of the phase diagram proposed in this paper and 
most of its attendent conclusions (with the speciflc exception of those which relate to the 
superconducting properties of this state) would remain valid . 

The multi-domain glass at intermediate H connects smoothly to the highly disordered 



glassy state obtained at large if, the putative vortex glass. The term "multi-domain" focuses 
attention on the intermediate level of translational correlations obtained within this phase 
in the interaction dominated regime. Melting on T scans appears to occur discontinuously 



out of the MG at these field values H3J4^]. Our MG phase is similar to the recently proposed 



"amorphous vortex glass" pB|J5(]| ]. Like the "phase-incoherent vortex glass" studied in detail 
in Ref. our MG phase is topologically disordered at the largest length scales [l51 



However, we emphasize that spatial correlations in the MG phase can be fairly long-ranged, 
unlike correlations in typical amorphous phases and in the disordered liquid; also our MG 
phase is a genuine phase, separated from the disordered liquid everywhere by a line of 
phase transitions. In addition, the phase diagram we propose differs significantly from ones 
proposed in earlier work. 

The multi-domain character of this phase is argued to manifest itself in the fact that the 
transition between BrG and DL phases can occur via several intermediate stages, correspond- 
ing to the melting of different domains as a consequence of disorder- induced inhomogeneities 
in the melting temperature This can result in large noise signals associated with 
a small number of fiuctuators PB| , ^ , intermediate structure in the ac susceptibility as a 
function of T f^JS^ , a stepwise expulsion of vortices |Q , strong thermal instabilities 



and a host of other phenomena specific to this phase |^5[, all features of the experimental 
data, which are hard to rationalize in other pictures. 

We argue that the MG phase melts via a first-order phase transition to the DL phase 
on T scans at intermediate H. Fig. 2 shows this first-order transition line terminating ||2^ 



at a tricritical point, where it meets a (continuous) glass transition line. (This topology 
is suggested by many experiments |jT^; however, recent experimental work, discussed in a 
later section, suggests a more complex topology, involving a critical end point, for the MG- 



DL transition line.) Some features of the phase diagrams proposed in Ref. ||3^ and in Fig. 
2, notably the two-step character of the transition out of the BrG phase in some regions 
of parameter space and the reentrant nature of the phase boundaries at low fields, were 



described in earlier work 31-34 



This paper addresses the following question: What links the phase behaviour of low and 
high-Tc systems? Physically, it is reasonable to expect that high-T^ and low-Tc supercon- 
ductors should have qualitatively similar phase behaviour. (This expectation motivates some 
of the features of Fig. 2.) However, Figs. 1 and 2 have some qualitatively different features: 
Fig. 2 implies that the Bragg glass always melts directly into a disordered glassy phase upon 
heating. This glassy phase is a continuation of the disordered phase obtained on cooling at 
high-field values. A subsequent transition, at intermediate H, separates the MG from the 
disordered liquid, as in Fig. 2. In contrast, consensus phase diagrams for high-Tc systems 
such as BSCCO and YBCO envisage at least one "multicritical point" (labelled (Hcr,Tcr) in 
Fig. 1) on the BrG boundary in the H — T plane. Below this point the BrG melts directly 



into the disordered hquid [|,P,|6,17,20 



This apparent distinction between the phase behaviour of low and high-T^ materials is 
surprising. This distinction is also manifest in the difference in phase behaviour exhibited by 
relatively pure samples of high-Tc materials which show an apparent multicritical point as 
discussed above, and somewhat more disordered samples. Even among the low-T^ systems, 
the data on relatively clean samples do not rule out a possible multicritical point, in contrast 



to the case for more disordered samples Mi 



Thus, at least two possibilities suggest themselves. We refer to these as scenarios (1) 
and (2). In the first scenario, a genuine distinction exists between more disordered and 
less disordered samples, leading to the following picture for the experiments: If the disorder 
is sufficiently weak the sliver of MG phase vanishes altogether and the BrG melts directly 



into the DL phase ||5^. If the existence of a low-field glassy phase is assumed (as shown in 
Fig. 2), there should be two multicritical points on the BrG-MG transition line, where the 
high-field and low-field glass transition lines meet the first-order BrG-DL phase boundary. 
As disorder is increased, these two points approach each other on the BrG-DL transition 
line. They then merge at a critical level of disorder. For stronger disorder, the BrG-MG 
transition line splits off from the MG-DL transition line, the two glassy phases are linked 
smoothly and the phase diagram of Fig. 2 is obtained. 



This scenario as well as the one outlined below assume the existence of a low-field glassy 
phase just above the lower critical field H^. Evidence for this phase is strongest in the 



low-Tc materials, in particular 2H-NbSe2, although there have been recent reports of 



similiar low-field glassiness in the high-Tc material YBCO |]58| , |59[| . The possibility of a low- 
field glassy phase has been raised theoretically |]BD|; simulations provide some support to 
this hypothesis ||6l]]. Recent experiments indicate that the domain of this low-field glassy 



phase can extend far above Hd in disordered samples of 2H-NbSe2 [04|,|62| . In addition, the 



transition between the low-field glassy phase and the BrG phase appears to be sharp ||62 |. 



The second scenario is the following: Any distinction between samples with different 
levels of disorder is only notional and the generic phase behaviour is that shown in Fig. 2 
|63| . Thus, the BrG phase would always be expected to transform first into a disordered 
MG phase and only then into the DL phase on T scans. It would then be necessary to 
rationalize the experimental observations of a possible single transition at intermediate H. 
A-priori, both scenarios are acceptable. However, they cannot both be correct if the phase 
diagram of a type-II superconductor with weak point pinning is universal. 

In this paper, it is suggested that the second scenario is a theoretically attractive and 
experimentally viable alternative to the first. This paper discusses links between the phase 
diagrams of Figs. 1 and 2 in the context of the available experimental data, simulation 
results and theoretical understanding. The phase diagram of Fig. 2 is proposed to be the 
more basic one; Fig. 1 is argued to be recovered from Fig. 2 when the BrG-MG and MG- 
DL transition lines shown in Fig. 2 come so close as to be separately unresolvable. (An 
alternative phase diagram which retains this feature but exhibits a more complex topology for 
the MG-DL transition line is discussed in a subsequent section.) It is argued that this is the 
case when thermal averaging over disorder is substantial, leading to an effective smoothening 
of the disorder potential. This discussion is motivated by the following conjecture for phase 
transitions out of the BrG phase: The BrG phase is conjectured to always transform directly 
into the MG phase and never in a single step into the DL phase, at any finite level of disorder 
|63| . Evidence from simulations, experiment and theory which supports this conjecture is 
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discussed. 

This paper contains six further sections. Section II summarizes the proposals of Ref. 
p5| for low-Tc systems and then supplements them with proposals appropriate for the de- 
scription of high-Tc materials. Quantitative estimates are then provided in support of these 
ideas. Sections III and IV describe a theoretical framework within which these proposals, 
particularly with regard to the phase boundaries, can be quantified - Section III discusses 
a Lindemann parameter-based approach to obtaining the phase boundary which separates 
BrG from MG phases while Section IV discusses a semi-analytic approach to the MG-DL 
transition based on a replicated liquid state and density functional theory. Section V com- 
pares the predictions of these sections with experimental data, principally on low-T^ systems. 
Section VI discusses the proposals of Section II in the context of the experimental data on 
both high and low-T^ systems. The concluding section (Section VII) links these results with 
previous work and suggests several experiments which could explore the ideas described 
here. 



II. UNIVERSAL PHASE DIAGRAM FOR 
DISORDERED TYPE-II SUPERCONDUCTORS 



The proposals of Ref. |^ derived from an analysis of peak effect systematics in the low- 



Tc materials 2H-NbSe2, CeRu2, Ca3Rh4Sni3, Yb3Rh4Sni3 and YNi2B2C. Ref. extending 
earlier work ||3^ , pointed out that these systematics suggested a general link to an underlying 
static phase diagram. The peak effect (PE) itself refers to the anomalous increase in jc seen 
close to Hc2{T); this increase terminates in a peak, as H (or T) is increased, before jc 
collapses rapidly in the vicinity of Hc2 . 



The results of Ref. relating to the low-Tc materials studied there are summarized in 
the first three proposals listed below; the reader is referred to the original paper for more 
details. These proposals are extended to high-Tc systems in this paper. The approach 
is primarily phenomenological - it draws on published experimental data, simulations and 
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theoretical input in attempting to describe a complete framework in which a large body of 
data on phase behaviour in the mixed phase with quenched point pinning disorder can be 
systematized and understood. 



A. Proposals 

1. The high-field "vortex glass" survives at intermediate field values as a sliver of dis- 
ordered, glassy phase intervening between quasi-ordered (BrG) and disordered liquid 
(DL) phases [^. (We use the phrase "multi-domain glass" (MG) to describe both 



the intermediate sliver phase and the putative vortex glass phase to which it connects 
smoothly.) Thus, in going from Bragg glass to liquid on increasing the temperature at 



intermediate field values, two phase boundaries are generically encountered [pl| - |33| , |67 
The first separates the BrG phase from the MG phase while the second separates the 
MG phase from the DL phase [pl| - |33| , |35[] . At low field values, the BrG-MG phase 



boundary is reentrant 0,^,^,0,0. For materials with intermediate to high levels 
of disorder, the glass obtained at high field values is smoothly connected to the one 
obtained at low fields ||35|. At low levels of disorder, the width of the sliver regime at 



intermediate fields is extremely narrow [pl| ,p4 



2. The sliver of glassy phase at intermediate field values defines the peak regime ||36|| , 
the regime seen between the onset of the increase of the critical current jc (at Tpi on 
temperature scans) and its peak value (at Tp on temperature scans) as functions of 
the applied field H and temperature T in materials which exhibit a sharp peak effect 
[0,0. The dynamical anomalies seen in experiments in this regime (slow dynamics, 
history dependence, metastability, switching phenomena etc.) arise as a consequence 



of the static complexity of the underlying glassy state P5[. Linking these anomalies 
to an underlying equilibrium transformation of the fiux-line system into a glassy state 
{i.e. the statics as opposed to the dynamics) suggests a simple and general solution of 
the problem of the origin of PE anomalies. 
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3. Structure in the sliver phase obtained at intermediate fields resembles a "multi- 
domain" structure, essentially an arrangement of locally crystalline domains The 
possibility of such an intermediate state is attractive because it represents a compro- 
mise between the relatively ordered lattice state and the fully disordered fiuid. At 
the structural level, it is a "Larkin domain solid" , originally believed to represent the 
fate of a crystal on the addition of quenched disorder ||69|. This proposal rationalizes 



a large body of data on PE anomalies including "fracturing" of the fiux-line array 
Q , the association of thermodynamic melting with the transition into the disordered 



liquid Eq] as well as conjectures regarding the dynamic coexistence of ordered and 



disordered phases in transport measurements in the PE regime [[70 |. 



4. The transverse size of a typical domain Rd in such a "multidomain" structure can 
be much larger than the mean-interparticle spacing a. It is conjectured that Rd is 
bounded above by Ra, the length-scale at which thermal and disorder induced fluctu- 
ations of the displacement fleld (calculated within the BrG phase) are of order a 



R.D is also conjectured to be largest in the intermediate-fleld, interaction-dominated 
regime but should decrease rapidly for much larger or smaller H, reflecting the in- 
creased importance of disorder both at high-field and low- field ends [^,35,|^. Similar 



statements should hold for the longitudinal sizes of such domains. In such a phase, 
translational correlations would be solid-like out to a distance of order R^j, but would 
decay rapidly thereafter. For sufficiently weak disorder, the size of a typical domain is 
large, refiecting the value of R^. In contrast, in the DL phase significant spatial corre- 
lations exist only only up to a few (1-2) intervortex spacings. Thus, the DL and MG 
phases at intermediate H should differ substantially at the level of local correlations 
for sufficiently weak disorder. 

5. At intermediate values of H, the glassy phase to which the Bragg glass melts is con- 
fined to a narrow sliver which broadens out both at much larger and much smaller 
H. Thermal smoothening of disorder reduces the width of this sliver. It can render 
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the sliver unobservable if the disorder is sufficiently weak, yielding an apparent single 
melting transition out of the ordered phase. This can happen for two reasons: First, 
domain sizes can be comparable to typical sample dimensions if the effective disorder 
is weak enough. Second, changes in critical currents signalling a transition between 
different phases can be too small to be resolvable; estimates also suggest that magne- 
tization discontinuities across the BrG-MG transition are very small. The apparent 
vanishing of the sliver explains the putative multicritical point |1T6| , pO| invoked in the 



context of phase diagrams for the mixed phase. While the sliver need not be resolvable 
in some classes of experiments which probe phase behaviour (for instance, it may not 
show up in dc magnetization measurements), it may be apparent in other experiments, 
particularly those which probe local correlations. For more disordered superconduc- 
tors, whether of the high-T^ or low-T^ variety, a two-step transition should generically 
be seen. 

6. A vortex glass phase of the type proposed by Fisher, Fisher and Huse may not exist. 
Gauge glass models for the VG transition do not show a finite T transition in three 
dimensions, once effects of screening are included |]7T|. Also, recent work questions 
earlier reports supporting a scaling description of the experimental data on disordered 
cuprates ||72,73]. Theoretically, the evidence appears to favour a low-temperature 



glassy phase without long-ranged phase correlations, at least in the weak disorder 



limit 1 11]. It is thus unlikely that the glassy phase obtained at low temperatures 



and elevated field values is a true phase- coherent vortex glass, in the sense of the 



original proposal [|r^. However, many experiments in the peak effect regime of low 
and high-Tc materials show that the multi-domain glass which intervenes between the 
BrG and the DL phases has strongly glassy properties PT|P^ , |5B|J^ , [T^ -[r^ . The data 



on low-Tc systems favours an equilibrium transition for the following reasons: (i) 
the approach to the intermediate regime is discontinuous and the phase boundaries 
reproducible, (ii) behaviour in this regime is dynamically very anomalous and, (iii) this 
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phase connects smoothly to the high-field (VG) phase, argued extensively to be glassy 
in nature. This suggests that the intermediate multi-domain glass is a true equilibrium 
glassy phase for reasons which may have nothing to do with the original vortex glass 
proposal. The precise nature of the symmetry breaking which distinguishes the MG 
phase from the DL phase is at present unclear; a "hexatic glass" of vortex lines and 
a "replica symmetry broken glassy phase" in the sense of modern approaches to the 
structural glass transition are both attractive possibilities. 

7. As a consequence of such glassiness, a hierarchy of long-lived metastable states can 
be obtained over significant parts of the phase diagram, leading to slow and non- 



trivial relaxation behaviour |[7^ , [77[| . Such states may dominate the experimentally 
observed behaviour at low temperatures and high fields and mimic the behaviour 
expected of an equilibrium vortex glass phase. As a further corollary, experiments at 
low temperature and high field values may often be probing strongly non-equilibrium 
behaviour, once time-scales for structural relaxation exceed experimental time-scales, 
thereby obscuring the underlying equilibrium phase transitions. 

B. Evidence for a Multi-domain state in the MG phase 

A substantial body of data favours a multi-domain structure in the regime intermediate 
between BrG and DL phases. A representative cross-section of this evidence is presented 
here. 

Recent neutron scattering measurements on disordered single crystals of Niobium 
(Hc2(4.2_ft') = 4.23 kOe,Tc{0) ^ 9K), which show a sharp peak effect, suggest fairly ordered 
crystal-like states obtained by heating following cooling in zero field, even for temperatures 
in the regime between Tpi (To in Ref. |78|) and Tp, the "peak regime" |7^. Following 



Ref. this regime defines the sliver of MG phase which intervenes between BrG and 
DL phases at intermediate H. In contrast, the metastable disordered states obtained by 
field-cooling through the peak regime show nearly isotropic rings of scattering, indicating 

15 



relatively short translational correlation lengths, comparable to those in a disordered liquid 
state. Experiments see a 30% reduction in rocking curve widths of Bragg spots between field 
cooled (FC) and zero-field cooled (ZFC) configurations demonstrating a fair degree of local 
order in the latter. That the relatively ordered states in the peak regime obtained on zero 
field cooling are actually lower in free energy can be demonstrated through experiments in 



which field-cooled states below Tp are annealed through the application of an ac field [ffS 
These data are consistent with SANS data on 2H-NbSe2 |^ . 



Recent muon-spin rotation experiments UTOipOj on BSCCO show a narrow intermediate 



phase where the asymmetry parameter associated with the field distribution function n{B), 

a = (Afi3)i/3/(A52^i/2^ (2.1) 

where AB = B — (B), jumps from a value of a ~ 1.2 at low T to a somewhat smaller value 
a ~ 1 before a further jump to a value close to zero at the irreversibility line. Large 
values of this asymmetry {a ~ 1.2) indicate a fairly well-ordered lattice structure, while 
values close to zero indicate a highly symmetric arrangement of vortices, as in a liquid. 
The observation of such asymmetric linewidths in this narrow intermediate phase can be 
associated to the multidomain structure of the sliver of MG phase expected to intervene 
between BrG and DL phases. These results would suggest a structure with a fair degree of 
local order, sufficiently different from the disordered liquid. 

As T is raised further, individual domains can melt due to disorder induced variations 
in the local melting temperature Tm, leading to a mosaic of solid-like and liquid-like (amor- 
phous) domains, similar to a regime of two-phase coexistence. This picture rationalizes 
several observations in the literature. Experiments on YBCO find that the magnetization 
discontinuities and first-order behaviour associated with the melting transition in the inter- 
action dominated regime at intermediate H can be associated with the MG-DL transition at 



higher fields, once the BrG-MG line splits off from the putative multicritical point ||81| , |82 
In addition, there appears to be a critical point on this first-order line, beyond which these 
discontinuities vanish. The idea that structure similiar to two-phase coexistence may ac- 
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count for some properties of the vortex array just below the melting transition has been 



raised in earlier work 183,8 



These features are easily explained within the picture outlined here. The melting tran- 
sition in this scenario should generically be associated with the MG-DL shown in Fig. 2, 
although the finite width of the sliver regime implies that some part of the discontinuity in 
density between fluid and solid phases can be absorbed across the width of the sliver If 



Rdi the characteristic domain size, is of order or smaller than the translational correlation 
length at freezing ^m-, sharp signals of melting are no longer expected [^]. On the other 
hand, ii Rd ^ at the melting transition, the behaviour should be equivalent to that 
of the pure system. Assuming Rd ^ initially, as Rd decreases on an increase in H 
(proposal (4) and the results of Refs. p4|,0) it approaches ^m] the point where Rd ~ 
defines the point where magnetization discontinuities should cease to be seen and thus the 
experimentally observed critical point. 

The possibility that the vortex glass phase obtained from the BrG phase on field scans 
might have a relatively low equilibrium density of dislocations was proposed recently by 



Kierfeld and Vinokur to rationalize some of the anomalies mentioned above UM: related 



ideas in the low-T^ context appeared in Ref . M^ . We argue here that such a picture is also 



applicable to the temperature-drwen transition out of the BrG phase, as is clear from Fig. 
2. 

Hypothesizing such a multi-domain in the intermediate sliver regime is consistent with 
the suggestion that a "fracturing transition" from a solid with quasi-long-range order into 
a intermediate state with correlation lengths of some tens of interparticle spacings is re- 
sponsible for the unusual open hysteresis loops obtained on thermal cycling within the peak 
effect regime in 2H-NbSe2 The idea that the regime which interpolates between the 

low-temperature ordered phase and the high-temperature disordered phase might have a 
relatively large degree of translational order has been used to rationalize the difference in 
behavior of the ac susceptibility in ZFC and FC scans in studies of the peak effect in several 
low-Tc materials [^]. 
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A fourth piece of evidence favouring a multi-domain or "fractured" state comes from 
voltage noise measurements, as pointed out in Ref. ||35|. Experiments see a substantial en- 
hancement in noise within the peak regime [^,^], a feature also seen in the ac susceptibility 
noise |3^. This noise is profoundly non-Gaussian in a regime where the PE is strongest; 
the number of independent fiuctuators contributing to this noise turns out to be small 



Experiments indicate a possible static origin for this noise, as opposed to a purely dynamic 
origin such as the interaction of fluctuating flow channels |5^. This supports a picture of 
domains of mesoscopic size in equilibrium with each domain an "independent fiuctuator" 
contributing to the noise. 



Recent simulation work on a pancake vortex model [gS] indicates that the PE may be 
associated primarily with a decoupling transition in the c-axis direction, in which c-axis 
correlations of vortices become short-ranged. This author and collaborators have argued 
recently that modelling the MG phase in a pancake vortex model in terms of a random 
stacking of perfectly crystalline layers should yield a fairly high-energy situation, with an 
energy which scales linearly with the area of each plane [^]. It is reasonable that some 
fraction of this cost can be relaxed by allowing dislocations to enter into each layer, leading 
to a situation far more like the "multidomain" structure proposed in this paper than a 



decoupled pancake-vortex state 



Other simulations, such as those of van Otterlo et. al. [p6| , |87| are consistent with an 
intermediate field MG phase with a correlation length intermediate between the solid and the 
disordered liquid. These simulations see primary peaks in the structure factor appropriate to 
the MG phase which transform into rings only across a second transition into the DL phase. 
The existence of these peaks indicates a fair degree of translational order, comparable to the 
sizes of the system studied, although the dynamics changes abruptly across the BrG-MG 
transition. Interestingly, van Otterlo et. al. conclude that they cannot rule out the existence 



of a sliver of MG phase always preempting a direct BrG-DL transition |35 



18 



C. Domain Sizes in the MG phase 



The scale for is itself set by a combination of the Larkin pinning length scale, at which 
disorder and thermal fluctuation induced displacements are of order ^, as well as the inter- 
line spacing a ||ll,p!I|. These can be derived from B{vj^,z). The transverse Larkin pinning 
length Rp is defined through B{Rp, 0) ^ A second length-scale, the longitudinal Larkin 
pinning length L^, is defined through 5(0, L^) ~ At transverse (longitudinal) length 
scales of order Rp (L^), the disorder-induced wandering of the displacement field equals the 
size of a pinning centre {^) at T = 0. A third important (transverse) length scale is the 
scale Ra at which disorder and thermal fluctuation induced displacements in the direction 
transverse to the field become of order the inter- line spacing a: B{Ra) — a?. 

At finite temperatures, replacing by max(,^^, < >) in the equations defining Lp 
and Rp above defines length scales L\ and Rc- Here < > is the square of a "Lindemann 
length"; these quantities are directly related to the critical current. In the Bragg glass, B{r) 
shows the following properties: For r^, z <C Rc-, correlations behave as in the Larkin 
"random force" model i.e. B{x) ~ x^"''. At length scales between Rc and Ra (the scale 
at which disorder- induced positional fluctuations become of order a), the random manifold 
regime, B{x) ~ x^'*'''". The exponent Crm ~ (4 — (i)/6 ~ 1/6. At still larger length scales, 
B{x) ~ log(x), yielding an asymptotic power law decay of translational correlations |l7|-|9|,P 



The Larkin length scales Rp and and the Larkin volume Vc ~ RjLc are estimated as 
lU = 2^^C66C44/n/^ and Rc = \pl^ cj^^ c^l j np . Using these, the standard weak-pinning 
expression for the critical current density jc follows [jl| : 



B' Wc, 

Rc and for a relatively clean 2H-NbSe2 single crystal are estimated in the following 
way [l|]. We anticipate that the following condition is satisfied: RciL\ > A, implying that 
non-dispersive {k = 0) values of the elastic constants can be used in our estimates. Then, 

Jc ~ I^J-^ jdp, (2.3) 
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~ -Re (2.4) 
a 

Assuming values of — 7K, X ~ 700A (if || c), ^ ~ 70A and a flux-line spacing 



a(= ^/2$o/V3fi) ~ 450A at S ~ IT. For 10"^ < jjjdp < 10-^ we then obtain: 

o r 150 (jVj,, = 10-6); 



(2.5) 

" ' 15 Uc/jdp = 10-^) 



and 

^6 r 240 (je/jrfp = 10-6); 

- ~ , (2.6) 

« [24 (je/j,p = 10-4) 

The domain size i?z) has been conjectured above (see proposal (4)) to be of order Ra in 

the intermediate H, interaction dominated regime. Ra is estimated via 

Ra^Rci^Y^^-, (2.7) 

where Rc is the (transverse) Larkin pinning length scale PJ9|JTl[]. At typical laboratory fields 
~ IT, we estimate 

(^)i/C™ ^ (1^)6 ^ 7 X 10^ (2.8) 

This leads to 

10^ < Ra/a < 106, (2.9) 

for a system with the parameter values described above. These are overestimates; allowing 
for the wave- vector dependence of elastic constants and prefactors neglected in these simple 
order-of-magnitude estimates should reduce these numbers considerably, possibly by factors 
of 10 or more. 

Assuming a multi-domain structure within the peak regime, a typical scale for each 
domain of Ra ~ lO'^a, a ~ 400A at if ~ IT and a (longitudinal) Larkin length comparable 
to the size of the sample in the c-axis direction, one can obtain an estimate for the number of 
"independent fluctuators" discussed above in the context of the noise measurements. For a 
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crystal of typical transverse area A ~ 1mm x 1mm ~ 10^^ the number of such fluctuators 
Nf is estimated as 

Nf = - 10^ - 10^ (2.10) 

numbers small enough to yield the strong non-Gaussian effects seen in the experiments, 
particularly if one assumes that not all domains are active contributors to the noise signal. 

Thermal fluctuations smear the effective disorder potential seen by a vortex line. If the 
thermally induced fluctuation of a vortex line about its equilibrium position exceeds ^, the 
effective pinning potential seen by a pinning site is strongly renormalized. An estimate of 



this effect at the level of a single pinned line ||60| yields 



Up ~ UpiO) exp[-c{T/Tdp)% (2.11) 

where Up(0) is the pinning potential associated with a single pinning site at T = 0, c is a 
numerical constant and Tdp is a characteristic temperature scale, the depinning temperature. 
The value of the depinning temperature T^p can be estimated from 

Tdpr^iUXceoh^Y^'. (2.12) 

Up is a measure of the pinning energy per unit length, eg is the line tension of a single vortex, 
c is a numerical constant of order unity and 7 is the mass anisotropy. Independent measures 
yield Tdp ~ 25 - AOK for YBCO |22|Q|81] . 



The ratio R = Tj°"'/Tj^*^'^ is then a quantitative measure of the relative importance of 
thermal fluctuations in low and high T^ materials. We use the following values: 



and 



NbSe2 : Up = lOK/A, ^ = 70 A, A = TOOA, 7 = 5 (2.13) 



YBCO : Up = lOK/A, ^ = 20A, A = 1400A, 7 = 50, (2.14) 



and obtain R ~ 10^, suggestive of the relative importance of this effect to the high-Tc 
materials vis a vis its irrelevance in low-T^ systems except very close to Hc2. A more 
accurate estimate of the magnitude of this effect comes from computing the ratio 
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''phigh 'Y'high 

R^, 10^ (2.15) 

m c 



where is the melting temperature of the pure flux-hne lattice and we have approximated 
Tm ~Tc in both cases. 

That thermal fluctuations lead to a substantial reduction in the effectiveness of disor- 
der close to the melting line in YBCO and BSCCO is clearly evident from experimental 
measurements of the irreversibility line in single crystals of these materials |l8l| , |82| , |9T| . The 
irreversibility line at intermediate fields {H ~ 3 — IT) in YBCO actually appears to lie well 
helow the melting line in weakly disordered samples Thus, thermal melting occurs in 



the almost totally reversible regime. Precision measurements reveal a weak residual 
irreversibility, with attendent jc values close to the melting transition of about QAA/om?, 
comparable to that obtained in the purest samples of 2H-NbSe2, but at much larger tem- 
peratures; Tm ~ SOi^ in YBCO at ~ 5T, compared to Tc's of about 7K for 2H-NbSe2. 



Similiar results have been obtained for BSCCO at fields less than about 0.05T [Q. In 



contrast, experiments on the low-T^, materials discussed in Refs. [p4| , |35| indicate an irre- 
versibility line which lies above both Tp and Tp^ lines, and subdivides the disordered liquid 
further into irreversible and irreversible liquid regimes. In these systems, the evidence for a 
two-step transition appears unambiguous. 

Such anomalously sma// values of jc (jc is proportional to the width of magnetic hysteresis 
loops and vanishes in the perfectly reversible regime), should lead to anomalously large 
values of Re- Using jc/jo ~ 10^'', and ^ and A values appropriate to YBCO, yields a 
characteristic domain size of order 1mm, clearly comparable to typical sample dimensions. 
Although these estimates are approximate ones, the physical intuition should be robust: 
characteristic domain sizes associated with an intermediate MG phase in high-Tc materials 
can be fairly large. Accordingly, in relatively pure samples of a small size, only a single 
sharp melting transition may be obtained, justifying the putative multicritical point and 
the direct transition into the liquid phase obtained in some measurements. 

The estimates above indicate that the length Ra can be far larger than typical correlation 
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lengths in the disordered fluid at the melting transition ~ 2 — 3a. SANS experiments on 
Nb indicate translational correlations of this order |7^. Thus sharp signals of melting can 
be expected at the second transition line intersected on a T scan out of the BrG phase; such 
signals should cease when Rd ~ C,m- 

The phase boundary between a multidomain solid and the disordered liquid can thus 
be argued on physical grounds to have a critical point. This is analogous to the gas-liquid 
transition in the following way: There is no symmetry distinction at the structural level be- 
tween a multi-domain state and a DL phase; the multi-domain solid has neither long-ranged 
crystalline order nor the power-law correlations of a Bragg glass. However, if a thermo- 
dynamic transition into a glassy phase with a symmetry different from that of the DL (or 
multi-domain state) exists, the associated phase boundary cannot terminate except at T = 
or on another phase transition line. Fig. 2 accounts for such an equilibrium glass transition. 
The topology of the MG-DL phase boundary is consistent with recent experimental work on 
single crystals of YBCO [^] and data on low-T^ systems p^ , |35[] . However, no arguments 
appear to rule out alternative locations for this line of glass transitions such as the one pro- 
posed in recent simulation work [BDI or more complex topologies, such as a critical end-point 
for the MG-DL transition line. 



III. BRG-MG PHASE BOUNDARY FROM A LINDEMANN PARAMETER 

APPROACH 

The field-driven BrG-MG transition occurs even at T = 0, where it is driven solely 
by changing the effective disorder. It represents a transition into a multidomain state at 
intermediate H on T scans. Such a state is substantively different from the equilibrium 
disordered fluid in terms of its local correlations and average density, provided the disorder 
is sufficiently weak. In contrast, the second transition, between MG and DL phases, is a 
close relative of thermal melting in the pure system. 

The simplest way to estimate phase boundaries such as those shown in Figs. 1 and 2 
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uses a Lindemann-parameter based approach. The treatment of the BrG-MG transition 
hne described here uses arguments similar to those of Giamarchi and Le Doussal (GLD) 
PJT7[] . It differs from their approach in the following way: GLD study this transition only at 
T = (where thermal fluctuations are absent) and at high temperatures where it is argued 
that disorder can be neglected and an expression for the melting line in the pure system 
used. The calculation here studies the crossovers in detail in the context of the experimental 
data on the low-Tc systems analysed in Ref. A more substantial difference between 



this approach and that of GLD is the proposal here that the high-temperature intermediate 
field transition out of the Bragg glass is not equivalent to thermal melting but rather to 
a transition to an intermediate 'multi-domain state" with the properties discussed in the 
previous section. 

A convenient phenomenological characterization of the melting transition in simple three- 
dimensional solids indicates that the transition occurs when the root mean-square fluctu- 
ation of an atom from its equilibrium position equals a given fraction (c^, the Lindemann 
parameter) of the interatomic spacing a. Giamarchi and Le Doussal suggested a possible 
generalization of this idea to the study of the instabilities of the Bragg glass phase 0. 
Taking over their proposal, we compute the ratio 

B{r^ = a)/a^ = c\, (3.1) 

(see Eq. |1.4|) , to obtain the melting line taking cj;, to be a universal number independent of 
H and T. 

In the random manifold regime, using a variational replica symmetry breaking ansatz, 
GLD derive the following relation 

B{r^) = ^h{r^/Ra), (3.2) 

where the crossover function 

h{x) ~ 1, X 1, (3.3) 

and 
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b(x) ~ x^O. (3.4) 

In a calculation which assumes wave-vector independent elastic constants, GLD obtain 
Rn as 



2a4 3/2 1/2 

Here C44 is the tilt modulus of the flux-line system and cqq its shear modulus, both calculated 
at zero wave vector. Up measures the pinning strength per unit length, a is the mean 
intervortex spacing and po is the areal density of vortices given by po — B/^o, with $0 the 
flux quantum: $0 = 2.07 x 10~^^Tm^. It is assumed that H ^ Hd, so that the effects of 
bulk distortions can be neglected. 

These expressions can be used to derive the full BrG-MG transition line in the H-T 
plane. For weak disorder, Ra ^ a. The following expressions for the elastic constants are 
used: 



C66 = ^(l-i^/5c2(r)r, (3.6) 



and 



- ^ B/B^T) ■ ^'-^^ 

Here eo = ($o/47rA^). The effects of a wave- vector dependent elasticity have been approx- 
imately accounted for in writing these expressions. Note that the elastic constants vanish 
as B ^ Bc2{T); at lower flelds the expression for Cqq should include a contribution from the 
line tension of isolated flux lines, which is neglected here. 
Some algebra then yields the intermediate expression 

where N is a. numerical factor. 

To simphfy U^/^^, we use the "core pinning" assumption 
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Up=P^e, (3.9) 

where p is a constant of order 1. This models the sources of pinning disorder as small-scale 
point defects on the scale of ^, which act to reduce Tc locally. Using 

H, = <^o/{2/{2)7:CX), (3.10) 

leads to 

I ^ 1 1 

^^'''^^(2^2M^/3^I^^'''^AV^' ^'-''^ 

Defining b = B/ B^ziT = 0) and assuming a temperature independent «;(= A/^), consid- 
erable algebra then yields the implicit expressions for the transition line h{t) given below, 
given the following assumptions about the temperature dependence of the penetration depth 
A and the upper critical field Hc2: 

1. Assuming a temperature dependence of A of the phenomenological two-fluid type i.e. 

X = X{T = {))/{l-t^fl\ (3.12) 

where t — T/Tc{H — 0) and a linear dependence of Bc2{T), i.e. 

B,2{T)^B,2{0){l-t), (3.13) 

yields the relation: 



(l-tl)2(l-t-6)7- 



(3.14) 



Here S is assumed to be constant; its value involves the product of the Lindemann 
parameter c^, and p in the following ratio: c^^k^ /p^, multiplied by a numerical 
constant of value 1.88 x 10^. Assuming k ~ 1 c^ ~ 0.22 and p ~ 1, yields E ~ 0.25. 
However, S is very sensitive to the value of cl used; changing by a factor of 2 to 
0.11 changes E by a factor of about 5 x 10"^. For cl — 0.15, k = 1 and p = 1, we 
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obtain S = 0.002. For this reason, we use S as a fitting parameter and consider values 
for it in tfie range 10"^ < S < 10^. 



Eq. ( |3.14| ) is plotted in Fig. 3, for different values of E, in the range 0.01 < S < 60; 



the MG phase lies above these lines while the BrG phase lies below them. Note that 
the transition lines are almost linear in (1-t) for large values of S. 

2. Assuming a temperature dependence of A of the phenomenological two-fluid type i.e. 

X = X{T = 0)/{l~t*y/\ (3.15) 

where t = T/Tc{H = 0) and a quadratic dependence of Bc2(T), i.e. 

5c2(T) = 5,2(0)(l-t2), (3.16) 

yields: 



2\6 



[i-t^)^{i-t^ -by 



(3.17) 



Again, the numerical value of S is dictated essentially by the value of c^, we choose a 
similiar range of S values as in the previous case. This relation is plotted in Fig. 4, for 
different values of S, in the range 0.01 < S < 50. These plots show a more substantial 
curvature than the analogous plots for case (1) above. 

In addition, assuming a temperature dependence of A of the following linear type A = 
A(T = 0)/(l — t)^/^ where t = T/Tc{H = 0) and a linear dependence of Bc2iT), i.e. 
Bc2{T) = -Bc2(0)(l — t) yields S" = ^j^^z^yr- This parametrization yield a BrG-MG phase 
boundary which is virtually a straight line, for the range of S values displayed in Figs. 3 
and 4. Note that all these results use different approximate parametrizations of Bc2{T)] 
the parametrizations involved in Eqs.( |3.1^ and ( p. 171 ) are both commonly found in the 



literature. Whether one or the other expression should be used should depend on the 
quality of fits to the experimentally obtained Bc2{T) line. We use the first parametrization 
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for the 2H-NbSe2 data (Section V), but the second parametrization for both the CeRu2 and 
the Ca3Rh4Sni3 data discussed in that section. 

For the phase boundaries to be physical, t + b must be less than or equal to 1. We perform 
the computations using Eq.( p.l4D . One point on the curve is (6 = 0, t = 1); 6 increases 
monotonically with 1 — t. At t = 0, the critical value of the magnetic field separating the 
BrG phase from the MG phase satisfies ^^^^^j^^Q-^y = S. For S small 6(0) ~ T,^^^, while for 
large S, 6(0) ~ 1 — The shape of the b — t phase diagram close to 6 = is also easily 

obtained. For t ^ 1, we have 6 ~ (1 — t). 

These calculations use a simple analytic parametrization of correlation functions in the 
BrG phase, obtained after many further approximations on an initially simplified Hamil- 
tonian, together with the restriction to a constant (field and temperature independent) 
Lindemann parameter. As a consequence, the quantitative predictions of this section should 
not be taken excessively seriously. However, the qualitative trends in the data appear to be 
borne out in this scheme of calculation, as discussed further in Section V. 

IV. SEMI- ANALYTIC CALCULATION OF THE MG-DL TRANSITION LINE 

This section discusses a calculational approach to the MG-DL phase boundary. The 
method outlined here uses results from a replica theory proposed by this author and Das- 
gupta (Ref. ||2^) for the correlations of a vortex liquid in the presence of random point 
pinning ||T^. This work examined the instability, within mean-field theory, of the DL phase 
to a crystalline state. 

The estimates obtained above indicated that typical domain sizes in the MG phase, 
could be much larger than the translational correlation length at freezing in systems with 
low levels of pinning. In addition, the transition out of the DL phase at intermediate H 
appears experimentally to be first order. A natural first approximation is to analyse, in 
mean-field theory, the instability of the liquid to an ordered crystalline state. This approach 
should thus represent the physics of the first-order part of the MG-DL transition line in the 



28 



interaction dominated regime at intermediate H. The issue of the nature of glassiness in this 
phase cannot be addressed by these methods. However, these calculations should provide a 
useful upper bound on the location of the actual MG-DL transition line. 

The calculations of Ref. ||2^ applied to a model for BSCCO which considered only the 
electromagnetic interactions between pancake vortices moving on different layers, ignoring 
their (far weaker) Josephson couplings |Q . Such a model is exact for a layered system in the 



limit of infinite anisotropy. For large but finite anisotropy, its predictions are quantitatively 
fairly accurate. It is argued here in some detail at the end of this section that the generic 
features of the results are relevant for the far more isotropic materials discussed in this 
paper. 

The analysis of Ref. was based on the replica method applied to a system of 



point particles interacting via the Hamiltonian 

H = Huneuc + ^ E ^(1 - I) + E M^i), (4.1) 

where V{r) is a two-body interaction potential between the particles and Vrf(r) is a quenched, 
random, one-body potential, drawn from a Gaussian distribution of zero mean and short 
ranged correlations: [V'ci(r)Vrf(r')] = K{\ r — r' |), with [■ ■ ■] denoting an average over the 
disorder. 

Using [InZ] = lim„_^o[(^" ~ l)/""-]? one obtains, prior to taking the n — > limit, a 
replicated and disorder averaged configurational partition function of the form 

„ n n N N 

, 7z- \ Ilc/rf exp( 

Here a, /? are replica indices and 1/°^(| - |) = - \)5a,p - I3K{\ r° - |) p7 



"1 "1 7T. Ti iv iV 

^/n<exp(-^EEEE^"^(|rr-r,^l))- (4-2) 



Eq. ( [4.2| ) resembles the partition function of a system of n "species" of particles, each 
labeled by an appropriate replica index and interacting via a two-body interaction which 
depends both on particle coordinates (rj,rj) and rephca indices (a,/3). This system of n 
species of particles can be treated in liquid state theory by considering it to be a n-component 
mixture and taking the n — limit in the Ornstein-Zernike equations governing the 
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properties of the mixture. These equations involve the pair correlation functions h"^ and 
the direct correlation functions C"^ of the replicated system. Assuming replica symmetry 
C"/3 = C^^)5o,p + - 5c,p) and = h^^^da^p + h^'^\l - 5o,p). 

Ref. obtained h^^\p,nd) and h^'^\p,nd) for the layered vortex system in BSCCO 
using the HNC closure approximation. To estimate the inter-replica interaction f3K{p,nd), 
the principal source of disorder was assumed to be atomic scale pinning centers such as 



oxygen defects ||90|. (Similar point defects are believed to act as the principal sources of 



disorder in the systems discussed in this paper.) A model calculation then yields 

(3V^'\p,nd) = -/3Kip,nd) ^ -F' exp(V/e')5n,o; (4.3) 

where pV^'^\p,nd) = (3V°'^{p,nd) with a ^ (3. ^ ~ IsA is the coherence length in the ah 
plane, and V ^ lO^^F^ for point pinning of strength dr^H^/STc, with d the interlayer spacing 
(~ ISA), r = pd^l/ATcX"^ and p = l/ksT. Defect densities of the order of W^/ciir' are 
assumed in the calculation of the prefactor . 



The calculated correlation functions C'^^\r) and C^'^\r) are then used as input into an 
appropriately generahzed (replicated) version of the density functional theory p6| , p!9 



in order to examine the effects of disorder on the freezing transition. Applying the replica 
treatment leads to the following functional in the n — limit: 
An 



dr 



p(r) In ^ - 6p{r) 



-\jdvj dv'[C^'\\ r - r' I) - C^'W r - r' |)][p(r) - p,][p(r') - p,] + . . . . (4.4) 

It is assumed that the density field is the same in all the replicas (p"(r) = p(r) for all a). 
The uniform liquid density is p^. This density functional resembles the density functional of 
a pure system with an effective direct correlation function given by 

C^//(| r - r' I) = C«(| r - r' I) - C^^\\ r - r' |), (4.5) 

In mean-field theory, the freezing transition of the pure system occurs when the density 
functional supports periodic solutions 
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P(i") = exp(iG ■ r), (4.6) 



G 



with a free energy lower than that of the uniform fluid |]96| , [l9| . Here G indexes the reciprocal 



lattice vectors of the crystal and pa represent the order parameters of the crystalline state. 
The properties of the freezing transition are controlled by C'^-^-^{q), the Fourier transform of 
C^f-f[r). In a one-order parameter approximation, only C^^^ {q) evaluated at the wave- vector 
q corresponding to the nearest neighbour spacing a i.e. at g = = 27r/a is required as 



input. The one-parameter density functional calculation of Ref. indicated that melting 
occured when the peC^^^^q = qm) attained a value of about 0.8 ||19|. This result agrees with 
the phenomenological Hansen- Verlet criterion for the freezing of a simple two-dimensional 
liquid Freezing occurs when the structure factor 

S{q) = l/{l-pC'^ff{q)), (4.7) 

evaluated at g.^ ^-c S{qm) attains a value of about 5. 

The following feature of the replicated density functional provides useful physical insight: 
Since C'^'^\qm) > 0, and C'^^\qm) is always reduced (although weakly) in the presence of dis- 
order, the equilibrium melting line is always suppressed by quenched disorder. The analysis 
of Ref. p3 showed that this suppression was very weak at low H but systematically increased 



at large if, indicating the increased importance of disorder at high fields. 

The semi-analytic approach to the suppression of the MG-DL phase boundary outlined 
here uses the following ideas: 

1. The diagonal direct correlation function C^^\q) is very weakly affected by disorder and 
can thus be approximated by its value in the absence of disorder. 

2. The off-diagonal direct correlation function C^'^\q) is a strong function of disorder and 
of H, the magnetic field 

3. C*^^-*(g) is well approximated at q = qm = 2% /a by its value at g = 0. 

Since the Hansen- Verlet condition is satisfied along the melting line, the following is true 
at melting: 
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p,C'ff = p,{C^^\qm) + C^^\qm)) ^ 0.8. (4.8) 

Now C(2) is a sharply decaying function in real space; its support is which for typical 
H <^ Hc2, is far less than a^. In Fourier space, therefore, its value at g = is close to its 
value at g = 0. We can approximate [p5|,3 

C^^\r) c:^ -f3V^^\r). (4.9) 

This scales with temperature as F^, implying that 

PeC^'\qn.) - ^. (4.10) 

Note that C^^\qm) increases as B is increased or as T is decreased, as is intuitively reason- 
able. 

We turn now to C«(g„). C(i)(g^) increases with a decrease in T; reducing T increases 
correlations. The variation in C'-^^ is expected to be smooth. We can therefore write 

C(i) {qm;T- AT, B) = {qm-, T, B) - q{B, T) AT, (4. 11) 

where q{B,T) (g < 0) is a smooth function of B and T close to the melting line. We 
are trying to find the effects to first order of adding C^'^\ so {B,T) can be replaced by 
{Brn, Tm) iu Eq. 4.8 above and C'^^^qm] B„i, Tm) by its value at freezing for the pure system: 
C'^^qm', Bm,Tm) — 0.8. The further approximation of neglecting the B and T dependence 
of g i.e. q{Bm,Tm) — q, with g a constant can also be made; at melting this dependence 
should be small provided a <^ A. 

An approximate expression for the suppression of the melting line from its value 
{Bm,Tm) = {BmiT),T) cau now be obtained. Using gAT^ = pBm/T^, 

Ar„ = ^ = = (4.12) 

with m = p/q approximately constant. Here AT^ is the shift in the melting temperature 
induced by the disorder. This relation predicts a larger suppression of the melting line 



at higher fields and lower temperatures, precisely as in the work of Ref. |^ and in the 
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experimental data. Given a parametrization of the pure melting line, Eq.( |4.12| ) can be used 
to estimate the effects of weak disorder on this line. 

This result can be combined with results from a calculation of the melting line in the 
pure system to obtain a simple analytic formula for the MG-DL phase boundary. At low 
fields, a Lindemann parameter-based calculation of this phase boundary yields 

B^{T)=C{T-Tc)\ (4.13) 



where Tc is the critical temperature and C is a constant Coupled with Eq.( [4.12|) above, 
this yields, in the presence of disorder 



B^%T) = C{T + ^"'^^^ - T^y. (4.14) 

This relation is plotted in Fig. 5 for different values of m together with the melting line 
in the absence of disorder i.e. with m = 0. The data shown in Fig. 5 use a Tc value of 
7K, as appropriate for 2H-NbSe2, a prefactor C of 1.15 kG (chosen purely for illustrative 
purposes) and values of Cm = 0, 0.2 and 0.5. Note that the suppression of the pure melting 
line by disorder is very weak if m is small. This suppression becomes progressively large as 
m is increased, or alternatively, at a lower temperature (larger H) for given m. These results 



i2 



are consistent, both qualitatively and quantitatively with the numerical results of Ref. [[27 



They are also consistent with the experimental results of Ref. ||3J] which find that fits of the 
Tp line in disordered samples of 2H-NbSe2 to a Lindemann-type expression are accurate at 
low field values but become increasingly inaccurate at larger fields. At large H the Tp line is 
systematically suppressed to lower T vis a vis. the fit, this suppression becoming apparent 
at lower fields for more disordered samples. 

These results were motivated using a parametrization of vortex lines in terms of inter- 
acting pancake vortices, a description valid for highly anisotropic layered superconductors 
in which the c-axis coherence length is far smaller than the interlayer spacing. How do they 
generalize to interacting linesl The calculation itself uses two quantities - the direct corre- 
lation function in the fluid phase for the pure system, evaluated at a wave-vector equal in 
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magnitude to the first reciprocal lattice vector of the crystal Gi i.e. C(| fc_L | = | Gi |, fc^ = 0). 
This quantity accounts for the tendency towards ordering at the length-scale of the inter- 
vortex spacing. The calculation which lead to Eg. ( [4. 121) above used only two properties of 



this correlation function: (i) the pure system freezes when the correlation function evalu- 
ated at this wave-vector achieves a particular (universal) value and (ii) the variation of this 
quantity upon a reduction of temperature can be parametrized simply, essentially linearly. 
Both these features should apply to the case of interacting lines. Numerical simulations of 
interacting lines provide some evidence for similiar quasi-universal attributes of the melting 



transition as seen in two and three dimensions, both in simulations and in experiments p9 

The other quantity required in the calculation, C'-^^ derives from properties of the cor- 
relation function of the disorder potential K{x — x'). This is independent of the (line or 
point) nature of the model for the fiux-line system used. Thus, it is apparent that the result 
obtained in Eqn. ([4.14|) , should be a fairly robust one which does not depend in detail on the 
fact that it was derived using results obtained from a model of interacting pancake vortices 
with only electromagnetic interactions. 

V. COMPARISON TO EXPERIMENTAL RESULTS 



Section III obtained expressions for the B — T phase boundary separating BrG and 
MG phases given Hc2(T) and A(T). In Section IV, qualitative and quantitative arguments 
were provided for the effects of weak disorder on the MG-DL transition line. This section 
compares the predictions of these sections with some of the available data on the phase 



boundaries in the experimental systems studied in Ref. ||3^ as well as data on the high-T^ 
material Nd2-:,Ce^CuOj, (NCCO). 

We do not attempt fits either to the YBCO or BSCCO data (although experiments 
relating to these are discussed in some detail in the following section), primarily because 
the results derived above ignore the effects of thermal renormalization of disorder, which are 
not negligible in these materials. (Effectively, these should lead to a strong T dependence of 
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S, which cannot be treated as a constant.) We have also not attempted to provide highly 
accurate fits to the experimental data in any of the cases below, since we merely wish to 
demonstrate that reasonable agreement with data can be obtained within the theoretical 
framework described in this paper. 

• 2H-NbSe2: Fig. 6 plots the experimental data for {Hpi,tpi) in single crystals of 
2H-NbSe2, taken from Fig. 2(a) of Ref. ||3^ together with a best fit line based on 



Eqn.( ^n^ ). As argued here (also see Ref. 0]), this line represents the BrG-MG phase 



boundary; tpi is Tpi/Tc.) An upper critical field of Hc2(0) = iGkG is assumed in the 
normalization of b. Assuming /t ~ 16, p ~ 1 and cl ~ 0.22, yields S ~ 60.0. The 
data are represented by an almost straight line over this field range, with deviations 
to the tune of about IkG appearing in the lower temperature range spanned by the 
data. The linear behaviour of Bc2(T) assumed in the use of Eq.( |3.14]) , is a feature of 
the experimental data in this field and temperature range. 

CeRu2: Data points representing the BrG-MG phase boundary, extracted from plots 
of {Hpi,Tpi), in Fig. 3(b) of Ref. are shown in Fig. 7, together with a best fit 



based on Eqn.(^]T^). An upper critical field of 68kG and a T^ of 6.2 K is assumed in 



the normalization to aid comparison to experimental data. Values of cl ~ 0.18, k ~ 1 
and p ~ 1 are assumed, yielding S ~ 0.01. 



Ca3Rh4Sni3: Points extracted from the plots of {Hpi,tpi) in Fig. 3(c) of Ref. ||34 



are shown in Fig. 8, together with best fits based on Eqn. (|3.17| . An upper critical 



field Hc2 of 45 kG is assumed in the normalization of the y— axis in the comparison 
to experimental data. Values of cl ~ 0.18, k ~ 1 and p ~ 1 are assumed, yielding 
S ~ 0.01. 



NCCO: Fig. 1 of Ref. [|123[j illustrates a very successful phenomenological fit to the 



experimental data on the BrG-MG phase boundary in NCCO. This phase boundary 
was assigned to the locus of onset points of the fishtail anomaly. We plot the line 
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obtained using Eq. p.l7| with S = 0.001, multiplied by an appropriate prefactor chosen 
to allow both curves to overlap as far as possible together with this fitting form = 
Bo(l — {T /TcYY^"^ in Fig. 9. (The somewhat smaller value of S here should reflect the 
increased role of disorder in these systems, as reflected in the parameter p.) Eq. |3.17 
predicts a marginally smoother and slower variation of the transition curve than the 
fitted form, but the qualitative trends appear to be the same. As a cautionary note, 
the identification of the BrG-MG phase boundary with the onset curves of the fishtail 
effect is not universally accepted. 

VI. DISCUSSION OF EXPERIMENTS 

To what extent are either of the two scenarios outlined in the Introduction supported by 
the experiments and the simulations? To support our proposal of the second scenario as a 
generic one for type-II superconductors with point pinning disorder, we take the following 
approach: First, we will show that the situation in which a sliver of intermediate phase 
intervenes between the BrG and the DL phases is actually far more common than earlier 
realized. We do this by showing that more recent experimental data on the high-T^ materials 
which initially showed a single melting transition strongly favours a two-step transition. 
Second, we will show that the first of the two transitions out of the BrG phase illustrated 
in Fig. 2, the "fracturing transition" into the MG phase, may not show up at all in many 
types of experiments commonly used to probe phase behaviour, such as measurements of 
the dc magnetization. 

This second point is illustrated through the following estimate: If a finite density of un- 
bound dislocations pd enters the sample at the BrG-MG phase boundary, the magnetization 
discontinuity across this boundary should scale, for small ~ 1/ i?^, as 

AM ~ AMoia/RaY, (6.1) 
where AMq is the magnetization jump in the pure system at the melting transition. Even 
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if Ra/a ~ 30, the corresponding induction jump AB ~ AM is of order 10 ^Mq or within 



noise levels in a typical experiment |15-5 



Some of the discussion will focus on the peak effect in critical currents seen close to the 



Hc2 phase boundary in weakly disordered systems. Within a simple Bean model for the 
critical state, the real part of the complex susceptibility, x', obeys 

X' ~ (6.2) 

once the ac field has penetrated fully within the sample. Here /? is a geometrical quantity 
related to the shape of the sample and h^c is the amplitude of the ac field. Thus, increases in 
the critical current density jc translate to reductions in x'. Transport measurements access 
jc more directly, although the presence of Joule heating in the transport measurements is 
often a significant factor, particularly in the peak regime where thermal instabilities are 
strong. 

The observation of an anomaly in magnetic hysteresis in the high-T^ superconducting 
oxides has attracted much attention in the past decade. This anomaly, an increase in the 
width of the magnetization hysteresis curve with field and a concomitant increase in jc 
(following an initial increase and subsequent decrease due to complete field penetration, the 
"first peak") is the "second peak" or "fishtail" anomaly [p^lOOj . The relation between jc 
and the width of the hysteresis loop follows from: 

jcoc AM = M{H 1) - M{H i), (6.3) 

the difference in values of the magnetization on the increasing and decreasing branches of 
the hysteresis loop. While several explanations have been proposed for this phenomenon, 
this behaviour is now believed to be correlated, at least approximately, with the field-driven 
transition between two vortex solid phases, identified by a large body of work as BrG and 
MG phases p2| . 



As pointed out in Ref. pH], the very structure of the phase diagram of Fig. 2 mandates 



the following: The peak effect in T scans at intermediate values of H should evolve smoothly 
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and continuously into behaviour characteristic of the BrG-MG field driven transition, both 
at low and high H . This is a simple consequence of the absence of a multicritical point. 
Insofar as the fishtail effect indicates a field-driven transition from a BrG phase to the MG 
phase, signatures of the fishtail phenomena should connect smoothly to signals of the peak 
effect on T scans. Thus, it is natural to argue for a connection between fishtail effects (in 
H scans) and peak effects (in susceptibility and transport measurements on T scans); both 
are phenomena which refer to the same underlying phase transition out of the Bragg glass 
phase. 

What is less clear is the precise connection of the fishtail feature to the underlying 
transition. There are at least three distinct characteristic fields associated with the fishtail 
feature. These are the onset field Hq (also called Hmm), the "kink" field (also called 



H3), the field value at which the magnetization shows a kink, and the peak field Hp ||103 
(also called H^p). The temperature dependences of these features can be quite different. For 
example, in YBCO, the onset is most often a monotonically decreasing function of T, while 
the kink and peak features can show non-monotonic behaviour. It has been argued that the 
kink feature most probably signals the underlying transition ||103|| , although there does not 
appear to be universal agreement upon this point. 

A. Low-Tc Systems 

This subsection summarizes features of the data on three low-T^ materials: 2H-NbSe2, 
CeRu2 and Ca3Rh4Sni3. Some aspects of this discussion have already appeared elsewhere 
35| , so the discussion here is brief. 



1. 2H-NbSe2 

Banerjee et. al. have reported data on ac susceptibility measurements on single crystals 
of 2H-NbSe2 of differing purity The purest crystal, X, has a Tc of 7.2K, while a crystal 



of intermediate purity, Y has a Tc value of about 7K. A more disordered crystal, Z, has a Tc 
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value of about 6K. The relative purity of the crystal was inferred from the values of jc/jdp 
for these materials. 

In sample X, Tp (where jc first begins to increase) and Tpi (where jc peaks) are almost 
coincident. The data indicate an extremely sharp peak, with a width comparable to or even 



smaller than the width of the superconducting phase transition in zero applied field P^ , |51 
The two-step nature of the transition is obvious in the sample with intermediate level of 
disorder, while this is a prominent feature of the most disordered crystal. The locus of Tp 
and Tpi lines, as functions of H, almost coincide in X, are resolvable as separate lines in Y 
and are clearly visible as separate features in Z. In sample Z, the Tp and Tp; lines clearly 
move apart at higher field values, suggesting that their ultimate fate would be to expand 
into the broad regime of vortex glass expected at high field values |P5|| . 

In the most disordered samples of 2H-NbSe2, onset and peak values of the peak effect are 
clearly and separately distinguishable down to low field values, where the Tp line begin to 
turn around, signalling the appearance of a reentrant glassy phase This feature 

suggests that the glass at high fields and the glass at low fields are smoothly connected, at 



least in samples with intermediate to high levels of disorder P3|. As Ref. p3 indicates. 



for purer samples, the extreme narrowness of the peak makes resolution of the features of 
the intermediate sliver region very difficult, although the peak itself is obtained to low field 
values. 

Ravikumar and collaborators see a discontinuity in the dc magnetization across the 
peak regime in single crystals of 2H-NbSe2 supporting the existence of a melting transition. 
The experiments could not resolve whether this feature was connected to Tp or to Tp; 
104|| . There is some evidence that the magnetization anomaly is connected to Tp from 



the combined magnetization and ac susceptibility measurements of Saalfrank et. al. ||105 
These experiments see a sharp peak in the derivative dM/dT, most likely associated with a 
step in the magnetization as a function of T at the location of the Tp line. Saalfrank et. al. 
note that the main peak of dm/dH showed no reentrant behaviour while the peak in the 
ac susceptibility shifted to lower temperatures at very low applied field values, supporting 
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the earlier proposal of a reentrant transition to a low field glassy phase |5^. As argued 
here, magnetization discontinuities signalling the transition to the liquid should generically 
be associated with the Tpi line. 

A substantial body of work exists on anomalies associated with the peak effect regime, 
particularly in the case of 2H-NbSe2 pT| , |5^ , |5B| , |55| , |T^ -[r5[] . Ref. ||^ proposed that these 
anomalies should be primarily associated with the static properties of the intervening vortex 
glass (equivalently, multi-domain) phase. This idea is discussed briefly here: (i) In the 
experiments, the transition to the anomalous peak effect regime on temperature scans from 
within the BrG phase occurs discontinuously via an apparent first-order phase transition. 
Thus the anomalies seen in the peak regime cannot be linked to the dynamics of plastic flow 
of a vortex lattice or even quasi-lattice state since the lattice phase is not linked smoothly to 
the intervening sliver phase (ii) the fact that peak-effect related anomalies in the depinning 
of a vortex lattice are only seen in the peak regime indicates strongly that a static phase 
transition may be involved and the solution cannot be traced to the dynamics alone. This 
point is stressed since many T = simulations of plastic flow (which do not have a phase 
transition) appear to be able to reproduce (some but not all) features of the peak effect 
phenomenon | p.06| ]. Finally, (iii) the concept of "softening" of elastic constants leading to a 
peak effect and associated anomalies may well be irrelevant; there appear to be few strong 
pre-transitional effects at the fracturing transition and the peak effect signal in x' jumps 
discontinuously in the experiments. 



2. Ca3Rh4Sni3 

The Ca3Rh4Sni3 system was first studied by Tomy et.al., who drew attention both to the 
prominent peak effect seen in this system and the related compound Yb3Rh4Sni3 as well as 
to the substantial similarities obtained with peak effects in other superconducting systems, 
notably UPd2Al3 and CeRu2 p9| , |i0[] . Tomy et. al. obtained H*, the field at which the onset 
of the peak occured both from susceptibility and dc magnetization experiments. Plots of H* 
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vs. T, show that the transition curve (which is interpreted here as signalhng the BrG-MG 
phase transition) appears to connect to the Hc2 hne only at T ^ T^. The region of fields 
and temperatures above H* was found to show strongly irreversible behaviour, in contrast 
to the behaviour below H*. 

Sarkar et. al. have studied in some detail the magnetic phase diagram of this compound 
through ac susceptibility measurements . Suceptibility data in this material also show the 



characteristic two-step feature exhibited by 2H-NbSe2 at intermediate fields. At low fields, 
the Tp and the Tpi lines come close but do not appear to merge at least until temperatures 
fairly close to Tc- The large field data are consistent with the bending backwards of the Tpi 
line indicating that the intermediate regime of vortex glass expands out at low temperatures 
and high fields, precisely as suggested in Fig. 2. 

3. CeRu2 

Early work on single crystals of the cubic Laves phase (C15) intermetallic superconductor 
CeRu2 by Huxley et. al. ||114| , p8| described a sharp reversible to irreversible transition in the 
mixed state close to B.c2 [|108|] . The possible relation of these results to a phase transition 
into a novel modulated superconducting phase (the FFLO state [|115|| ) was discussed by 



these authors and later investigated extensively by others. Careful longitudinal and Hall 
resistivity measurements by Sato et. al. have revealed that the peak effect survives in this 
material to far larger fields than the original measurements indicated ||117|| . Fig. 2 of Sato 
et. al. summarizes the results of their measurements and is to be compared to Fig. 2 of 
this paper; the Hp(T) curve in this figure is to be related to our Tpi{H). Note that this hne 
appears to extrapolate smoothly to Tc, although signatures of magnetic hysteresis are not 
observed below about 2T. Sato et. al. point out that their observations rule out the FFLO 
state as a possible explanation of the PE anomaly. These experiments suggest that pinning 
is substantially weakened at intermediate H, although signatures of a peak effect remain. 
Transport measurements on CeRu2 by Dilley et. al. show a striking peak effect signal. 
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together with a strong dependence of measured quantities on the magnitude of the transport 



current | |118|| . Dilley et. al. note the existence of strong thermal instabihties in the interme- 
diate peak effect regime, together with profound hysteresis. Many of these features are seen 
in 2H-NbSe2 |§5|. Muon-spin rotation studies of CeRu2 ||119|| have been invoked to suggest 



that a FFLO mechanism may be involved in this material. These measurements see a rapid 
decrease in the second moment < Ai?^ > of the field distribution function across a line 
in the phase diagram; such a decrease naively suggests a sharp increase in the penetration 
depth A. Yamashita et. al. point out that pinning induced line distortions should generically 
lead to broadening of the line-width and not to its reduction. However, a sharp transition 
to a multi-domain phase, with short-ranged correlations in the field direction would give 
essentially similar results to those of Yamashita et. al. [pT| . 

In a recent illuminating study of the peak effect anomaly in single crystals of CeRu2, 
Tenya et. al. have shown that the peak effect lies in the irreversible region of the phase 
diagram, counter to some previous work ||107|| . Their results establish the following: The 
onset of the peak, at Tpi, corresponds to an abrupt transition from a defective lattice to 
a regular one. The state of the sample in the peak regime is argued to be intermediate 
between an ordered structure and a fully amorphous one. These ideas are closely related to 
those presented here. 

Work by Banerjee et. al. shows a remarkable similarity between the vortex phase diagram 
in a weakly pinned single crystal of CeRu2 with that of more disordered samples of 2H-NbSe2 
||33| , |3^ . This close similarity has been used to argue that the origin of peak effect anomalies 
in these two materials may well be common [^,0 . 



B. High-Tc Systems 

This subsection summarizes data relating to peak and fishtail effects in four high-Tc 
materials: BSCCO, YBCO, NCCO and (K,Ba)Bi03. 
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1. YBCO 



Early torsional oscillator experiments of d'Anna et. al. measured complex response in 



the PE regime of untwinned YBCO | 120 |. A prescient comment of d'Anna et. al. regarding 
their data is particularly noteworthy, for its overlap with the proposals made here. d'Anna 
et. al. commented that ..the melting phenomenon, giving rise to the peak effect, the loss 
peak and the kink in resistivity in so-called clean YBCO, is in fact a complex, two-stage 
phenomenon.. |T^]. The first stage was ascribed to a lattice softening leading to a change 
in the pinning regime, as in conventional approaches to the peak effect, while the second was 
associated with a collapse of the shear modulus at melting and/or a decoupling transition. 
In transport measurements in fields upto lOT, d'Anna et. al. found two boundaries. 



between which noise due to vortex motion rose to a maximum ||113|| . The lower boundary 
was identified with the peak effect; the noise falls below the limit of resolution at another 
boundary, somewhat above the first one. These results are to be compared to those obtained 
in similiar transport and susceptibility measurements on 2H-NbSe2, which see a 
sharp increase in noise amplitudes only in the peak effect regime, with a discontinuous 
onset. 

Ishida et. al. have simultaneously measured ac susceptibility and magnetization and see 
a sharp peak effect in the field range 0.1- 1.5T | |121| |. The peak of the peak effect signal in ac 
susceptibility was found to correlate exactly to the location of the magnetization jump which 
signals melting. Ishida et. al. draw attention to a small dip feature in x' at a temperature 
Ts, preceding the peak and the associated magnetization discontinuity. This subtle feature is 
attributed to a "synchronization effect between lattice and pinning sites" , as a consequence 
of lattice softening. (The possibility of such synchronization has often been discussed in the 
past as a mechanism for the peak effect ||6^.) The close similarity of this phenomenology 



with that discussed for low-Tc systems in Ref. is emphasized here. Ishida et. al. also 
comment on the enlarged regime in which x' shows non-trivial signatures of the transition 
in comparison to the magnetization jump which occurs at a sharply defined temperature. 
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These results are very simply understood using Fig. 2. We would argue here for the 
natural identification = Tpi. The width of the transition region is related to the width of 
the sliver phase, argued here to be always finite in scenario (2), while the most prominent 
signatures of melting should generically be obtained across a line, the MG-DL transition 
line, which we suggest is the true remnant of the vortex lattice melting line in the pure 
system. 

Shi et. al. have observed a "giant" peak effect in ultrapure crystals of YBCO [ [122|] at 



fields ranging from 0-4T ||122|1 . This feature is most prominent at intermediate H, repro- 



ducing the phenomenology of the peak effect in 2H-NbSe2. Shi et. al. demonstrate that 
the jump in the magnetization coincides, to within experimental accuracy, with Tp. Shi et. 
al. also comment that peak effect signals can be seen in fields up to 6T, although only by 
going to far higher ac drive frequencies (1 MHz). This is consistent both with the weakening 
of the effective disorder in the intermediate field regime and the existence of a sliver phase 
with finite width. 

Rydh, Andersson and Rapp | |127| ], in transport studies of untwinned YBCO, obtain three 



regimes of flow: (1) a low temperature creep regime (2) an intermediate flow/ creep regime 
(3) a small regime in which the resistivity drops as a consequence of the peak effect and (4) 
a high-temperature fluid regime. Rydh et. al. find that the dip in resistivity is correlated to 
the melting transition. It is suggested that the width of region III reflects a distribution of 
melting temperatures in the solid and argue that the onset of melting occurs at the minimum 
of the resistive dip. In our interpretation, regime III would be assigned to the multi-domain 
regime. 



Suggestive work by Rassau et. al. [|8J] studies magnetization relaxation in the vortex 
solid phase of YBCO. Rassau et. al. argue that a well-defined region exists below T^ which 
can be quantified as a coexistence of solid and liquid phases and demonstrate that the vortex 
solid state can be pinned with different strengths, depending on prior history. The processes 
which lead to these different states arise across a narrow temperature region immediately 
below Tm. This is precisely the phenomenology indicated by experiments on 2H-NbSe2 and 
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related low-Tc compounds which show a peak effect; we draw the readers attention to Ref. 
1)771 in this regard. We would argue that the intermediate region interpreted by Rassau et. 
al. as relating to a regime of two-phase coexistence is our multidomain phase. 

Detailed insights into the vortex phase diagram of untwinned YBCO comes from the work 
of Nishizaki et. al., who present magnetic measurements on the vortex lattice in a clean 



sample |^2[- large H and low T, the data show a fishtail feature in the magnetization 
which narrows, as T is increased, to a bubble-like feature. Nishizaki et. al. discuss the 
existence of an anomalous reentrant behaviour in jc, inferred from their data via a Bean 
model-based calculation, for temperatures in the range 68/^ < T < 7AK. They find that jc 
drops very sharply at low fields to almost unobservable levels, then picks up to form a peak; 
while hysteresis decreases and disappears at T = 70K in the intermediate field region (3.5 
T< fioH < 6T), it reappears again and is maximum around 10. 5T. 

Nishizaki et. al comment that this reentrant magnetization is seen only in high quality 
samples with a lower pinning force and comment that these results indicate that the pinning 
force for the untwinned samples is remarkably reduced in the intermediate- field regime. This 
comment accords with our earlier discussion of the strong effects of thermal renormaliza- 
tion of disorder in the intermediate field, interaction dominated regime. The shape of the 
irreversibility line in the data of Nishizaki et. al. is particularly noteworthy in this respect. 

Very accurate ac susceptibility measurements using a local Hall probe have been per- 
formed by Billon et. al. in the field and temperature regimes in which Nishizaki et. 



al. find reversible behaviour and a single melting transition. These experiments find that 
the critical current is actually finite below the melting temperature but is extremely small 
(~ 0.4A/cm^), leading to an irreversible magnetization unresolvable by standard supercon- 
ducting quantum interference device magnetometry. These experiments see a peak effect 
below the apparent melting transition, a feature inaccessible in conventional SQUID based 
measurements. These results clearly illustrate that signals of a two-step transition in weakly 
disordered single crystals of high-T^ materials may be very hard to access, particularly if 
discontinuities in jc or magnetization across the first transition are small. 
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Nishizaki et. al. argue that in the high temperature region above Tcp, the BrG-MG 
transition hne (H*(T) in the notation of Nishizaki et. al.) turns down and meets the Hm(T) 
mehing hne helow Hep for irradiated YBCO. Thus they argue that the second peak effect just 
below Tm may be closely related to enhanced vortex pinning due to vortex lattice softening 
i.e. the conventional peak effect. The connection of this observation to the discussion of 



the phenomenology of the peak effect presented in Ref. is stressed here, as is the link 
outlined earlier between fishtail features and PE features, given the phase diagram of Fig. 
2. 

Nishizaki and collaborators have also reported results on the phase diagram of untwinned 
YBCO crystals on varying the Oxygen stoichiometry. For fully oxidized YBCO crystals 
(YBa2Cu30j^, y ^ 7, Tc ~ 87.5K), Nishizaki et. al. find that the first order melting transition 
can be tracked to high fields (upto SOT). For optimally doped YBCO (YBa2Cu30j^, y ~ 6.92, 
Tc — 93K) and slightly underdoped YBCO, it is suggested that a vortex slush phase can 
exist between the second-order and first-order transition lines. This "slush" phase intervenes 



between the glass and the liquid phases of the vortex system |^ . 

How is the slush phase to be understood using the ideas presented here? Fig. 2 presents 
a particularly simple topology for that part of the phase diagram which involves the MG-DL 
transition line, in which the first-order MG-DL transition becomes continous for sufficiently 
large H . Such a phase diagram does not contain a slush phase. However, the situation could 
be more complex. For example, the continuous part of the MG-DL transition line obtained 
at large H could meet the first-order part at a critical end point (CEP). The first-order 
transition above the CEP would then be interpreted as a liquid-liquid transition, between 
phases of different densities, with a further, symmetry- breaking transition into a glassy phase 
occurring at lower T . The phase with a larger average density would then represent the 



"slush" phase [109 



Abulafia et. al. ||110|| study magnetization relaxation in YBCO crystals, finding two 



distinct regimes of relaxation below and above the peak in the fishtail magnetization. These 
regimes were interpreted by these authors as signalling plastic vortex creep in the vortex 
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lattice state. This interpretation has been questioned by Klein et. al. ||111|| (see also Ref. 



112|| ) who argue that the fishtail effect marks a transition from a low-field ordered phase 
to a high-field glassy structure, an interpretation in agreement with ours. It is argued that 
topological order (in the glassy phase) can he quenched at large length scales, with diverging 



harriers as j goes to zero | |1 1 1|| , an idea in agreement with our suggestion that the multi- 
domain glass has fixed topological order unlike the liquid, but glassy attributes. With this 
interpretation of the data of Abulafia et. al. their phase diagram (Fig. 4 of Ref. [|11CI|] ), 
is very similar to the one we have proposed. In particular, the melting line and the glass 
transition line (if it is identified with the locus of fishtail peak positions), converge only at 

T ^ Tr. 



Recent experiments by Pissas et. al. [|126|| on single crystal YBCO samples probe some- 



what more disordered samples than those studied by Nishizaki et. al.. An abrupt change 
in slope of the increasing part of the virgin magnetization curves is seen. This feature is 
sharpest at intermediate T; it broadens or vanishes at low T and is hard to discern for 
T > 75K. This abrupt change of slope, associated with a curve H3(T), merges with the sec- 
ond peak feature B.sp{T) at a temperature T*. Both the irreversibility line and the melting 
line lie ahove the second peak line. The structure of this phase diagram is to be compared 
to the ones shown in Ref. for low-Tc materials in which the irreversibility line lies above 
the BrG-MG transition line. The phase behaviour depicted in Fig. 4 of Pissas et. al. agrees 
well with the ideas presented here as well as the structure of Fig. 2 with the identification 
of H3 as the BrG-MG transition line. 

To summarize, experiments on untwinned YBCO support the presence of an intermediate 
phase which intervenes in equilibrium between BrG and DL phases, particularly at low 
fields {H < 2 — AT). For larger H, a substantially reversible regime is entered, close to 
the putative melting transition. In this regime, precision experiments see weak residual 
irreversibility as well as a peak effect just below the melting transition. (SQUID-based 
magnetization experiments see only a single transition here.) It is natural to interpret these 
features as signalling the continuation of the sliver between high and low-field ends as well 
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as the profound weakening of disorder effects in the intermediate H, interaction dominated 
regime, as a consequence of the averaging of disorder by thermal fluctuations. We argue 
that this weakening of effective disorder is responsible for the apparent single-step transition 
experiments see at intermediate H. At still higher fields, the fishtail feature splits off from 
the melting line, yielding an expanding regime of vortex glass phase. The magnetization 
discontinuities associated with this melting line survive to very high fields (~ SOT) in the 
purest samples. For more disordered samples, the separation between the BrG-MG and the 
MG-DL transition lines is plainly apparent. 



2. BSCCO 

The earliest data to provide thermodynamic evidence of a sharp melting transition in 
single crystals of a high-Tc superconductor were the Hall probe measurements of Zeldov 
and collaborators on BSCCO ||T5|. These experiments obtained a discontinuity in the mag- 
netic induction sensed by Hall probes at the surface of a BSCCO platelet. Practically no 
hysteresis was seen at any one probe, while different probes showed slightly different values 
for the melting field, indicating a spread of melting temperatures within the sample. The 
temperature driven melting transition in relatively pure single crystals of BSCCO occurs in 
a highly reversible regime. As reported for YBCO, the irreversibility line actually lies within 
the domain of the solid phase for clean crystals. These experiments see a single melting 
transition for weakly disordered systems. These and related experiments supported earlier 
structural evidence, from neutron scattering |]14| and muon-spin rotation for a sharp 



transition out of a quasi-lattice phase on increasing both B and T. 

Very recent magneto-optic measurements by Soibel et. al. have reexamined this issue 
||5^ . These experiments see a global rounding of the transition as a consequence of quenched 
disorder, argued by these authors to be due to a hroad distribution of local melting tempera- 



tures at scales down to the mesoscopic scale . These experiments reveal a remarkable and 



complex coexistence of fluid and solid domains in the larger crystals with the properties that 
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while the local transition is sharp, the global solid-liquid transition is rounded by quenched 
disorder. Soibel et. al. found that the interfacial tension between solid and fluid phases 
was very small, indicating that the vortex melting process was governed to a large extent 
by disorder. 

Recent muon-spin rotation experiments, on three sets of crystals whose properties range 
from overdoped to underdoped, show very clear indications of a two step transition, as 
pointed out by the authors. The flrst transition on increasing T was associated with a 
decoupling transition between the layers. In contrast, in our picture, this would be the 
transition into a multi-domain phase. This identiflcation is supported by the observation of 
a values close to unity in this phase indicating a substantial degree of local order. 

Soibel et. al. see substantial "supercooling" across the flrst-order melting transition, with 
relatively small domains of fluid coexisting with domains of solid. The putative supercooled 
state then transforms abruptly into the crystal. We point out here that precisely such a 
scenario would be obtained in the context of the phase diagram of Fig. 2 with the single 
proviso that the intermediate "coexistence" regime which Soibel et. al. see would be assigned 
to our proposed "multidomain" state, an equilibrium phase with glassy properties. This is 
entirely consistent with our interpretation of the muon-spin rotation results of Blasius et. al.. 
We point out that from related studies on 2H-NbSe2, it is known that the intermediate state 
on fleld cooling is relatively highly disordered, with short-ranged correlations resembling 
those in the liquid. By analogy, one would expect similar disordered states to be seen in the 



BSCCO case on fleld cooling, precisely as seen |52[| . 

Ooi et. al. in a study of vortex avalanches in BSCCO through local magnetization and 
local permeability measurements, see a stepwise expulsion of vortices in a distinct temper- 
ature regime helow the flrst-order transition on decreasing temperature scans [Q. This 
expulsion ceases abruptly below a second, lower temperature called by these authors as T^^. 
In between Tm and Td-, the experiments see broad-band noise with a power law spectrum. 
This phenomenology precisely reproduces related noise data in the peak effect regime of 
2H-NbSe2 and related systems. 
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Ooi et. al. also report a possible single step transition for weakly disordered samples 
and a temperature- dependent peak effect, which they associate with inhomogeneities in the 
melting field, for more disordered samples ||129|| . These results, presented as schematics as 
Fig. 5 of Ooi et. al. have interesting parallels with the phase diagrams presented here. We 
would argue that the regime that Ooi et. al. assign to inhomogeneities in the melting field 
should be attributed instead to a genuine thermodynamic phase, the sliver of disordered 
MG phase, which intervenes between the BrG and DL phases. 

A simulation by Sugano et. al. finds clear evidence for a two-stage melting of the 
disordered vortex-line lattice in BSCCO ||130|| . These authors find that the low- field melting 
transition out of the Bragg glass occurs always via an intermediate "soft" glass phase; their 
phase diagram (Fig. 3 of Ref. ||130|| ) is to be compared to the one shown in Fig. 2. Their 
proposal that the glassy phase is subdivided into strongly pinned and weakly pinned regimes 
has interesting overlaps with our suggestion of the increased importance of non-equilibrium 
effects at low T. 



3. NCCO 

The onset of the fishtail peak in magnetization measurements on a Nd1.85Ceo.15 Cu04_5 
crystal, a layered cuprate system with substantial anisotropy and a relatively low T^ (~ 23 
K) has been probed via local magnetization measurements using an array of Hall probes 
||123|| . These experiments (in particular the phase diagram of Fig. 1 of Ref. [p.23|| ), show 
an onset field for the peak effect Ho which is always distinct from the irreversibility line 
separating an "entangled solid" phase from the vortex liquid phase. This phase diagram 
supports the ideas presented here, in particular the contention that a single sharp melting 
transition out of the BrG phase is not generic (as scenario (1) would suggest) and provides 
strong support for scenario (2). 
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4. (K,Ba)Bi03 



Small angle neutron scattering studies of the isotropic high-Tc material (K,Ba)Bi03 
(Tc ~ 30K) suggest a phase diagram which is very close to the one displayed in Fig. 2 of this 
paper |p.28|| . In this system, as in the NCCO system discussed previously, the transition line 
between the quasilattice and the glassy state always lies below the transition line separating 
the glass from the liquid. Klein and collaborators comment that the SANS data in this 
system support the connection of the fishtail effect with the transition to a glass out of the 
quasi-lattice BrG phase ||128||. In addition, another interesting feature of the data is the 



relatively smooth decrease in the intensity of a Bragg peak on field scans, while the width of 
the peak remains virtually constant, suggesting a relatively high degree of local correlations 
in the putative vortex glass phase, as we would expect for a multi-domain solid. 

VII. CONCLUSIONS 

This paper has presented arguments in favour of a generic phase diagram for type-II 
superconductors with quenched point pinning disorder. This phase diagram. Fig. 2, differs 
from others proposed earlier. We suggest that the relatively ordered BrG phase generically 
transforms into an intermediate glassy state with solid-like correlations out to relatively large 
length scales rather than directly into a liquid. We have pointed out that many experiments 
are, in fact, consistent with this proposal and discussed how signals of the first transition, 
from the BrG to the MG phase, may often be hard to detect. We have also presented a 
simple analytic parametrization for the BrG-MG phase boundary, drawing on earlier work 
as well as provided an analytical expression for the MG-DL phase boundary. These are 
consistent with the experimental data. 

Hypothesizing a multi-domain state in the intermediate- if regime is consistent with the 
experimental data. This suggestion rationalizes the association of thermodynamic melting 
with the MG-DL transition line; for a prior suggestion in this regard see Ref . |2^ . In addition. 
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given the expectation that the effects of disorder increase as H increases, it indicates a 
physical mechanism for the occurence of a critical point on the first-order MG-DL transition 
line. 

Related theoretical and simulation work which bears on ideas proposed here include 
work by Feinberg | |125| |, Ikeda [ [I24| ] and van Otterlo et. al.. Ikeda has suggested that a 
vortex glass instability may accompany a first-order solidification transition in the presence 
of weak point pinning disorder [p.24|] . Feinberg has argued that melting of the Bragg glass 



may be to an intermediate glassy phase, via a reentrance of single particle pinning ||125 



van Otterlo et. al. point out that their simulation work cannot rule out the possibility of a 
sliver of glassy phase always intervening between ordered and fluid phases, as in the phase 
diagram of Fig. 2. Several authors have commented [ |128|| , on the relation of these simulation 



results to their experimental data. Some of the ideas presented here also have overlaps with 



recent work by Kierfeld and Vinokur |2J]. A recent comprehensive survey of the status of 
vortex glass phases |Tl| makes much the same points as we do regarding the absence of a 
true phase- coherent MG phase as envisaged in the original proposal 0. Several of the ideas 
presented here draw from extensive work on peak effect phenomena in 2H-NbSe2 and related 
materials, summarized in Ref. [Q. 

Many of the proposals presented here are, in fact, experimentally testable. One is our 
proposal of a multi-domain phase. Experiments which probe local correlations should be 
able to access such structure. We have suggested that in the regime where the sliver of 
intervening glassy phase narrows, structural correlations in this intermediate phase should 
become large. This proposal is testable both in simulations and in experiments. 

The second is our identification of the first of the two phase transitions out of the BrG 
phase on increasing T as a thermodynamic phase transition. Signals of this phase transition 
in the form of entropy jumps or singularities in the specific heat will be extremely small, 
since they reflect ordering at the scale of units of 10^ vortex lines or more. However, it 
remains to be seen whether high-precision experiments might be able to resolve our sug- 
gestion of two equilibrium, thermodynamic phase transitions generically separating the low 
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temperature Bragg glass phase from the disordered hquid. Such experiments would provide 
crucial evidence in favour of the ideas presented here. 

Perhaps the single most novel proposal of this paper is the conjecture that the low- 
T, quasi-lattice Bragg glass phase in all superconductors should generically melt into an 
intermediate glassy phase before finally transforming into a liquid This possibility 

violates none of what is known about the phenomenology of the experimental data, nor the 
simulations nor available theoretical evidence. Further tests of the ideas presented here as 
well as first-principles calculations of the phase diagram of Fig. 2 would be very welcome. 
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FIGURES 




TEMPERATURE (K) 



FIG. 1. The current view of the phase diagram of type-II superconductors ||£.17,16,ll,pO,24|, 
incorporating the effects of quenched point pinning disorder and thermal fluctuations. In addition 
to Meissner and normal phases, this phase diagram subdivides the mixed phase into three phases 
- the Bragg glass (BrG), the vortex glass (VG) and the disordered liquid (DL). The BrG melts 
directly into the DL phase through a first-order melting transition on T scans at intermediate H. 
On H scans at fixed low T, the BrG phase transforms discontinuously 128,29) into the disordered 



VG phase. The VG phase is understood to transform via a continous phase transition, with the 
exponents linked by relations calculated in Ref. |^, into the DL phase. This line of continuous 
phase transitions meets the BrG phase boundary at a "multicritical" point, labelled {Hcr,Tcr) 
in the figure. The first-order direct transition from BrG to DL may persist beyond the putative 
(RcrjTcr) |4^,B2,2^] terminating in a critical point (not shown in this figure, see Refs. [p^ , |TT| ). 



The line which separates the DL from the normal phase is a crossover line and not a true phase 
transition. This phase diagram does not show the small regime of reentrant glassy phase expected 
at low field values; for a phase diagram which includes this see Ref. [p^. 
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FIG. 2. Proposed universal phase diagram for type-II superconductors incorporating the 
effects of thermal fluctuations and quenched random point pinning. The term multi-domain glass 
(MG) is used here in preference to "vortex glass" ; for a discussion of this point and of the properties 
of this phase, see text. The other phases are described in the caption of Fig. 1. Note that the MG 
phase intrudes between BrG and DL phases everywhere in the phase diagram. At intermediate 
H where interactions dominate, the MG phase is confined to a slim sliver but broadens out again 
for sufficiently low H. The putative "multicritical point" {}lcr,Tcr) of Fig. 1 is identified with 
the location where the BrG-MG and MG-DL transition line first approach so closely that they 
cannot be separately resolved. The inset to the figure expands the boxed region shown in the main 
panel at low fields and temperature values close to Tc(0) and illustrates the reentrant nature of the 
BrG-MG phase boundary at low fields. The transition line separating the MG from the DL phase 
is first order at low fields but terminates at a tricritical point, where it meets a line of (equilibrium) 
glass transitions. The transition out of the BrG phase is generically first-order. 
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FIG. 3. Plot of Eq. 3.14 , the expression for the BrG-MG phase boundary in the b,t plane (6 
and t are the reduced field and temperature H/Hc2 and T/Tc respectively), derived in Sec. Ill, 
for different values of S, in the range 0.01 < S < 60 (see text). E is a phenomenological fitting 
parameter involving the Lindemann parameter and fundamental constants. 
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FIG. 4. Plot of Eq. 3.17, the expression for the BrG-MG phase boundary in the b,t plane 



and t are the reduced field and temperature H/Hc2 and T/Tc respectively), derived in Sec. Ill, 
for different values of S, in the range 0.01 < S < 60 (see text). S is a phenomenological fitting 
parameter involving the Lindemann parameter and fundamental constants. 
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FIG. 5. Plot of Eq. 4.14, the expression for the MG-DL phase boundary incorporating the 
effects of quenched disorder as described in the text, derived in Sec. IV. The data use a Tc value of 
7K, as appropriate for 2H-NbSe2, a prefactor C of 1.15 kG (chosen purely for illustrative purposes) 
and values of Cm = 0,0.3 and 0.5. Note that the pure melting line is increasing suppressed 
by quenched disorder, this suppression becoming larger as the applied field H is increased, in 



agreement with the predictions of Ref. |27| and experiments. 
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FIG. 6. Plot of the experimental data for {Hpi,tpi) in single crystals of 2H-NbSe2 (points taken 
from Fig. 2(a) of Ref. [Q) together with a best fit line based on Eqn.(^rT4). Here tpi is Tpi/Tc- 
An upper critical field of Hc2(0) = 4:6kG is assumed in the normalization of b and S = 60 is used; 
see text. 
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FIG. 7. Data points representing the BrG-MG phase boundary in CeRu2 (points taken from 
plots of {Hpi,Tpi) in Fig. 3(b) of Ref. together with a best fit based on Eqn.( ^rT7| ). An upper 

critical field of 68kG and a Tc of 6.2 K is assumed in the normalization of the y— axis. A value of 

S = 0.01 is used; see text. 
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FIG. 8. Data points obtained from plots of Hpi,tpi in Ca3Rh4Sni3 (from Fig. 3(c) of Ref. 
|) are shown together with best fits based on Eqn.( ^!T4 . An upper critical field Hc2 of 45 kG is 
assumed in the normalization of the y— axis in the comparison to experimental data. A value of 
S = 0.01 is used; see text. 
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FIG. 9. The solid line is the solution of Eq.|3T7| with S = 0.001. The dotted line is the fitting 



form Bm = Bq{1 — {T /Tc)^)^^"^ , which provides an accurate fit to the onset value of the fishtail 
effect in NCCO []123| . The prefactors are chosen to allow the two curves to overlap as extensively 
as possible. 
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